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The one-dimensional spin-| XY model was introduced in 1961 by Lieb, Schultz and Mattis. 1 
The authors of that celebrated paper found that several statistical mechanics calculations for 
that spin model could be performed exactly because it can be rewritten as a model of non- 
interacting spinless fermions with the help of the Jordan- Wigner transformation. Evidently, 
the formulation of the spin-| XY chain in terms of fermions permits one to give a magnetic 
interpretation to the results derived for one-dimensional tight-binding spinless fermions. As an 
example of where this relationship has been exploited one can cite the papers on the Lloyd 
model 2 ' 3 and corresponding papers on the spin-| XX chain in a random Lorentzian trans- 
verse field. 4 ' 5 The work reported in the present paper has been inspired by some results on 
one-dimensional tight-binding Hamiltonians for periodically modulated lattices 6 ' 7 and spinless 
Falicov-Kimball model. 8 ' 9 Combining the approach developed in those papers and the treat- 
ment of the Lloyd model presented in Refs. 2, 3 we shall calculate exactly the random-averaged 
one-fermion Green functions (that yield the density of states and therefore the thermodynam- 
ics) for the periodic nonuniform spin-| XX chain in a random Lorentzian transverse field. We 
shall treat a few particular chains in order to discuss the changes in the density of states and 
magnetic properties induced by periodic nonuniformity and diagonal disorder. It should be 
noted, in passing, that the periodic nonuniform spin- 1 ; XX chain was considered in several 
papers dealing with the spin-Peierls instability in a spin-^ XX chain 10,11 (see also Refs. 12- 
16). However, those papers concentrated on the influence of the structural degrees of freedom 
upon the magnetic ones, rather than on the properties of a magnetic chain with regularly al- 
ternating exchange couplings. One should also mention a study of a spin-| XX model on a 
one-dimensional superlattice 17 (such a model can be viewed as a nonuniform chain with period- 
ically varying exchange coupling) but consideration was restricted to the excitation spectrum. 
Our communication is also related to the work in Ref. 18 and may be viewed as a further 
study of the effects of periodic nonuniformity and randomness on magnetic properties of spin-| 
chains. 

Let us consider a cyclic nonuniform XX chain of N — > oo spins \ in a transverse field 
described by the Hamiltonian 

N N 

h = E nx + 2 ]T / n + 

n=l n=l 
N , 1\ N 

= ^2^n[ S+S- - - J + ^ I n (s+SrT+l + S n ■ (!) 

71=1 V Z/ 71=1 

Here Q n is the transverse field at the site n, and is assumed to be a random variable with the 
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Lorentzian probability distribution 

Qon is the mean value of the transverse field at the site n, and T n is the width of its distribution. 
2I n is the exchange coupling between the sites n and n + 1. After making use of the Jordan- 
Wigner transformation the model is recasted into a chain of spinless fermions governed by the 
Hamiltonian 

N , i \ N 



H 



( 1 \ 

2 fin ( C+C n - - J + T n (4c n+1 - C n C+ +1 ) (3) 



n=l x 2/ n=1 



(the boundary term that is non-essential for the calculation of thermodynamic quantities has 
been omitted). Note, that for the non-random case (T n = 0) assuming the transverse field in (3) 
to be uniform one arises at the Hamiltonian considered in Ref. 6. In addition, after substitution 
f2 n — > Uw n , I n — > —t, Eq. (3) transforms into the Hamiltonian of a one-dimensional spinless 
Falicov-Kimball model in the notations used in Refs. 8, 9. In model (3), considered here, not 
only the transverse fields that are independent random Lorentzian variables having the mean 
values and widths of their distribution, but also the exchange couplings between neighbouring 
spins, vary from site to site. 

Let us introduce the temperature double-time Green functions 
Gf im (t) = ^(±t)({c n (t),c+(0)}), Gl m (t) = (l/27r)r oo d^exp(-i^)GJ m (^±i e ), e - +0, 
where the angular brackets denote the thermodynamic average. Utilising a set of equations for 
GJ im = G^ m (u ± ie), and performing random averaging using contour integrals 2-5 , one finds 
the following set of equations for the random-averaged Green functions 



[U ± ir n — Q 0n ) G^ m — I n -lG n -i jTn — InGn + l, m ~ $nm- (4) 



Here (...)= Jf^ dfiip(Qi) . . . Jf^ dfijvp(fiAr)(- • •)• Our task is to evaluate the diagonal random- 



averaged Green functions G^ n , the imaginary parts of which give the random- averaged density 



of states p(u), 



1 N 



PM = ^I>G* (5) 



that in its turn, yields thermodynamic properties of spin model (1). It is a simple matter to 



obtain from Eq. (4) the following representation for Gf m 

1 



G± 



"" u±\T n -Sl Qn - A--A+' 



T 2 

1 n-l 



CO ± ir„_i — "0,n-l 



I 2 

1 n-2 



A+ = 

n 



uj ± ir„_i_i — Q 



ra+1 — "0,n+l 



/ 2 



(6) 



LJ±ir n + 2 — f2(),ri+2 ~ 



Eqs. (5), (6) are extremely useful for examining thermodynamic properties of periodic nonuni- 
form spin-i XX chain in a random Lorentzian transverse field when periodic continued fractions 
emerge. 

Consider at first, a non-random case. For a regular alternating chain 

fliIiQ 2 I 2 Q 3 I 3 QiIifl 2 I 2 Q 3 I 3 • • • one generates periodic continued fractions having a period 3, 
and after some calculations one gets 



p(co) 



0. 



1 \X{w) 



3tt 



/cm' 



if to < c 6 , c 5 < to < c 4 , c 3 < to < c 2 , Ci < a;, 
if c 6 < to < c 5 , c 4 < to < c 3 , c 2 < < ci; 



#(u;) = if + J 2 2 + J 3 2 -(co- Qi) (w - Q 2 ) - (w - fii) (w - ft 3 ) -(to- Q 2 ) (to - Q 3 ) , 



C(co) 



AT 2 T 2 T 2 



Il (to - Q 3 ) + Il (to - no + 7 3 2 (w - Q 2 ) 



fiO (W - fi 2 ) (W - ft 3 ) 
6 

= - n - co 



(7) 



Here Ci > . . . > c% denote six roots of the equation C(to) = that can be found by solving 
two cubic equations. For a regular alternating chain S7 1 / 1 r2 2 -^2^3^3^4-^4^i-^i^2^2^3-^3^4-^4 • • • 
periodic continued fractions having period 4 emerge, and after some calculations one finds 



p(to) 



0, 



if to < d$, d 7 < to < d 6 , d$ < to < d 4 , d 3 < to < d 2 , d\ < to, 
if d 8 < to < d 7 , d 6 < to < d 5 , d A < to < d 3 , d 2 < to < d\\ 

w(to) = il(2to -n 3 - n 4 ) + i\(2to n 4 ) + ii(2u - n ± - n 2 ) + il(ito -n 2 - n 3 ) 



i 1WMJ 

4^^M' 



-(u 


— Vli)(uo — Q 2 )(to — fi 3 ) — (to — Qi)(to - 


-n 2 )(to 


- n 4 ) 


-(to- 


- — Jl 3 )(o; — fi 4 ) — (to — Q 2 )(to — 


n 3 )(to- 


-n 4 ), 


V(to) 


= Ailiplil - [(to - no - n 2 ) (to - 


n 3 ) (to 


-n 4 ) 




-il (to - n 3 ) (to - n 4 ) - il (to - 


^0 


-n 4 ) 




-i 2 (co - no (oo - n 2 ) - ^ 2 - 


n 2 ) (a; 


-n 3 ) 




+/ 1 2 / 3 2 + / 2 2 / 4 2 ] 2 = - 


-n(- 


-d,-)- (8) 
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Here, d 1 > . . . > d 8 are the eight roots of the equation V{uj) = that can be found by solving 
two equations of 4th order. Let us note that all Cj and dj are real since they can be viewed as 
eigenvalues of symmetric matrices. 6 

The simplest periodic nonuniform spin-^ XX chain in a random Lorentzian transverse field 
r2oiri/ifio2r2^2^oiri^i^o2r2-^2 ■ • • requires a calculation of periodic continued fractions with 
period 2. The random-averaged density of states for this case becomes 

i | y(d) | 



2tt B(u) 



\ /-j-, „ , B(u)+B'(u) nlu s /n „ „ s B(u)-B'(u) 



B'(u) 



[u - - q 02 ) - r 1 r 2 - if - jf] 2 - [(u - q 01 )t 2 + ^r!] 2 - 
B"{u) = 2\(u- n 0l )(co - n 02 ) - T l v 2 - 1 2 - 1 2 ] [(uj - n 01 )v 2 + n 02 )v l ] . (9) 



In principal, there are no difficulties in considering more complicated periodic nonuniform chains 
apart from the fact that the calculations become somewhat cumbersome. 

Let us turn to a discussion of the results obtained for the density of states. First, note 
that in the limit of a uniform transverse field and exchange coupling, Eqs. (7) and (8) re- 
produce the well-known result for the uniform spin-| XX chain in a transverse field: p{uj) = 
l/n^AI 2 — {u — VL) 2 if 4/ 2 — {uj — Q) 2 > and p{uj) = 0, otherwise. This result for the uniform 
chain also follows from (9) if Q i = ^02 = fi, Ti = T 2 = 0, I± — I 2 — I. The density of states 
(9) contains the result for uniform spin-| XX chain in a random Lorentzian transverse field 4 



p(u) = T(l/vr)Im [l/yj(u ±iT- fi ) 2 - 4/ 2 J if Q 01 = tt 02 = Q , T 1 = T 2 = T, h = I 2 = I. In 
addition, the density of states (9), in the non-random limit Ti = T 2 = 0, coincides with the 
density of states (8) with Q 1 = Q 3 , Q 2 = fi 4 , I± = I 3 , I 2 = J 4 , i.e. for a regular alternarting 
chain VliIiVl 2 I 2 VLiIiVL 2 I 2 ... as should be expected. 

In Figs. 1 - 3 we display the density of states, together with both the dependencies of the 
transverse magnetization 

1 roo ft 

^='2 loo dEp(E) tanh 2kf (10) 
on the transverse field at zero temperature, and the static transverse linear susceptibility 

on temperature, for few particular chains considered. Initially, let us discuss a non-random 
case. The main result of introducing the nonuniformity is a splitting of the initial magnon 
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band with the edges Q — 2 | 7 | , + 2 | 7 | into several subbands. The edges of the subbands 
are determined by the roots of equations C{uj) = 0, V{u) = 0. p(u>) is positive inside the 
subbands, tends to infinity inversely proportionally to the square root of u when u approaches 
the subbands edges, and is equal to zero outside the subbands. For special values of parameters, 
the roots of the equations that determine the subbands edges may become multiple, the zeros 
in denominator and numerator in the expression for the density of states cancel each other, and 
as a result one observes a smaller number of subbands. In Figs, la, 2a we show p(u) for some 
periodic nonuniform chains having periods 3 and 4 to demonstrate the energy band scheme 
in the presence of nonuniformity. Note that the splitting caused by periodic nonuniformity is 
not surprising, since the periodic nonuniform chain can be viewed as a uniform chain with a 
crystalline unit cell containing several sites of the initial lattice. On the other hand, one expects 
several subbands for a crystal having several atoms per unit cell (see for example Ref. 19). We 
now turn to the density of states given by (9). First, note that in the non-random case one 
finds two magnon subbands the edges of which are given by 

1 



{61,62, 63,64} 



n 1 + n 2 ± - n 2 ) + 4 (I h I ± I h \Y 



(12) 



.2 

(Fig. 3a). Introduction of the uniform diagonal Lorentzian disorder Ti = T 2 = T results in 
smearing out of the edges of subbands (Fig. 3b). For extremely nonuniform diagonal Lorentzian 
disorder r x = 0, T 2 7^ and with a small strength of disorder only one subband may be mainly 
smoothed, but with increasing of the strength of disorder both subbands become smeared out 
(Fig. 3c). 

The described dependencies of the density of magnon states on introduction of periodic 
nonuniformity and disorder affects the behaviour of thermodynamic quantities of the spin 
model. For example, the temperature dependence of the specific heat 

c = Jf^ dEp(E) (gjlr) / cos h 2 2§r ^ or a non - ran dom periodic chain may exhibit a two-peak 
structure consisting of low-temperature and high-temperature peaks. Let us comment in some 
detail on the magnetic properties of the spin chains considered. The splitting of the magnon 
band into subbands caused by nonuniformity has interesting consequences for those properties. 
Consider for example the transverse magnetization (10). Since for T — > 0, tanh ^ tends either 
to —1 if E < 0, or to 1 if E > 0, one immediately finds because of the appearance of subbands 
that the low-temperature dependence of m z against Q in the non-random case must be com- 
posed of sharply increasing parts (when, with increasing Q, E = moves inside each subband 
from its top to its bottom) separated by horizontal parts (when, with increasing of Q, E = 
moves inside the gaps). Evidently a number of plateaus in the low-temperature dependence of 
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the transverse magnetization on transverse field is determined by a number of subbands. Every 
cusp in the the dependence of the transverse magnetization on Q induces a singularity in the 
dependence of the static transverse linear susceptibility (11) on Q. One can easily show that 
the dependence against Q, at T — 0, is the same as the dependence of p(uo) against VL — uo. 
The latter dependence can be derived from the densities of states depicted in Figs. 1-3. The 
value of Yzz at T = in the temperature dependence of xll at f2 = is determined by the 
value of p(0). Therefore, the nonuniformity and randomness may essentially effect this value 
as well as the temperature dependence of the static transverse linear susceptibility. In different 
temperature regions one finds both an enhancement (and even a divergence, as in Fig. 3g), or 
suppressing of the curves — xlz versus T. This can be seen nicely in Figs, lc, 2c, 3g, 3h, 3i. 

To summarize, using continued fractions we have obtained rigorously the density of magnon 
states for periodic nonuniform spin-| XX chain in a random Lorentzian transverse field. Con- 
tinued fraction representation of the solution of Eq. (4) is extremely useful for calculation of 
the density of magnon states. The attractive features of this approach can be seen even for the 
uniform case. In this case one omits performing twice the Fourier transformation while solving 
Eq. (4) in a standard manner and evaluates straightforwardly the desired p(u). The advan- 
tages of the continued fraction approach become clear while treating periodic chains already 
with the smallest period of 2. Periodic nonuniformity leads to a splitting of the magnon band 
into subbands that in its turn leads to the appearance of new cusps and singularities in the 
low-temperature dependences of the transverse magnetization and the static transverse linear 
susceptibility on transverse field, respectively. In the random case the spectacular changes in 
these dependences are smoothed. Periodic nonuniformity and randomness may either enhance 
or suppress the temperature dependence of — Xzz- Changing the degree of periodic nonuni- 
formity one may to some extent influence the detailed shapes or Xzz against f2, or Xzz 
against T. The described approach may be of considerable use for examining simple models of 
spin-Peierls instabilities, especially in the presence of disorder. Analysis of the thermodynamic 
properties of non-random and random Lorentzian periodic nonuniform spin-| XX chain in a 
transverse field and its stability with respect to a lattice distorsion will be given in a separate 
paper. 

The author is grateful to O. Zaburannyi for discussions and J. W. Tucker for critical reading 
of the manuscript and useful comments. He acknowledges the financial support of the Research 
Assistance Foundation (L'viv). He is also indebted to Mr. Joseph Kocowsky for continuous 
financial support. 



7 



References 



[i] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7 
[8 

[9 
[10 

[11 
[12 

[13 

[14 

[15 

[16 

[17 

[18 



E. Lieb, T. Schultz, and D. Mattis, Ann. Phys. (N.Y.) 16, 407 (1961). 
P. Lloyd, J. Phys. C 2, 1717 (1969). 

W. John and J. Schreiber, Phys. Status Solidi B 66, 193 (1974). 

H. Nishimori, Phys. Lett. A 100, 239 (1984). 

O. Derzhko and J. Richter, Phys. Rev. B 55, 14298 (1997). 

S. W. Lovesey, J. Phys. C 21, 2805 (1988). 

Ch. J. Lantwin and B. Stewart, J. Phys. A 24, 699 (1991). 

J. K. Freericks and L. M. Falicov, Phys. Rev. B 41, 2163 (1990). 

R. Lyzwa, Physica A 192, 231 (1993). 

P. Pincus, Solid State Commun. 9, 1971 (1971). 

R. A. T. Lima and C. Tsallis, Phys. Rev. B 27, 6896 (1983). 

K. Okamoto and K. Yasumura, J. Phys. Soc. Jpn. 59, 993 (1990). 

K. Okamoto, J. Phys. Soc. Jpn. 59, 4286 (1990). 

K. Okamoto, Solid State Commun. 83, 1039 (1992). 

Y. Saika and K. Okamoto, |cond- mat / 9 5 1 1 14 . 



A. Fujii, |cond-mat/9707137 . 



L. L. Gongalves and J. P. de Lima, J. Magn. Magn. Mater. 140-144, 1606 (1995). 

O. Derzhko and O. Zaburannyi, J. Phys. Stud. (L'viv) 2, 128 (1998); 

O. Derzhko, O. Zaburannyi, and J. W. Tucker, J. Magn. Magn. Mater. 186, 188 (199£ 



[19] A. S. Davydov, Tjeorija tvjerdogo tjela (Nauka, Moskwa, 1976) (in Russian). 



8 




Figure 1: 
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Figure 2: 
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Figure 3: 
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List of figure captions 



FIG. 1. Density of states (a), transverse magnetization versus transverse field Q at T — (b), 
and static transverse linear susceptibility versus temperature at Q = (c), for the nonuniform 
chain QhQhnWIi^Wh ■ ■ • , h = 0.5. I 2 = 0.5, I 3 = 0.5 (solid lines); I 2 = 0.5, I 3 = 1 
(dashed lines); I 2 = 0.25, J 3 = 1 (dotted lines). 

FIG. 2. The same as in Fig. 1 for the nonuniform chain niiVlI 2 ni 3 Ql4QIiVlI 2 VlI 3 £ll4 ■ ■ ■ , 
i! = 0.5. I 2 = 0.5, I 3 = 0.5, h = 0.5 (solid lines); I 2 = 0.5, J 3 = 0.5, h = 1 (dashed lines); 
I 2 = 0.5, J 3 = 0.25, h = 1 (dotted lines); I 2 = 1, I 3 = 0.25, I 4 = 1 (dash-dotted lines). 

FIG. 3. Density of states (a,b,c), transverse magnetization versus transverse field Q at T = 
(d,e,f), and static transverse linear susceptibility versus temperature at Q = (g,h,i), for the 
nonuniform random chain £701X1/1^02^-^2^01X1/1^02^^2 ■ ■ ■ , h = 0.5, Oqi — ^01 — 0.5 + ^2. 
a, d, g: non-random case T 1 = T 2 = 0, I 2 = 0.25 (solid lines), I 2 = 0.5 (dashed lines), I 2 — 1 
(dotted lines); b, e, h: uniform disorder J 2 = 0.25, Ti = T 2 = (solid lines), Ti = T 2 = 0.5 
(dashed lines), T 1 = T 2 = 1 (dotted lines); c, f, i: nonuniform disorder I 2 = 0.25, Ti = 0, 
T 2 = (solid lines), T 2 = 0.5 (dashed lines), T 2 = 1 (dotted lines). 
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